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^D ! A simple exact analytical solution of the relativistic Duffin-Kemmer-Petiau equation within the 
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framework of the asymptotic iteration method is presented. Exact bound state energy eigenvalues 
and corresponding eigenfunctions are determined for the relativistic harmonic oscillator as well 
as the Coulomb potentials. As a non-trivial example, the anharmonic oscillator is solved and 



^ ■ the energy eigenvalues are obtained within the perturbation theory using the asymptotic iteration 

o. 

^ ■ method. 
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I. INTRODUCTION 

Exact analytical solutions to relativistic wave equations are important in relativistic quan- 
tum mechanics since the wave function contains all the necessary information to describe 
a quantum system fully. There are only a few potentials for which the relativistic Dirac, 
Klein-Gordon and Duffin-Kemmer-Petiau (DKP) equations can be solved analytically. So 

II nn 

far, many methods such as the super-symmetric (SUSY) p|, shape invariance [2, |3!], factor- 
ization and path integral |j, |5|, la |7[ etc have been developed to solve the relativistic wave 



equations exactly, or quasi-exact ly, for potentials like Coulomb, harmonic oscillator, Posch 
Teller and exponential type ones. In recent years, an asymptotic iteration method for solv- 
ing second order homogeneous linear differential equations has been proposed by Ciftci et 



al 
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10|. This method has been applied to solve the non- relativistic radial Schrodinger 
and Dirac equations for various potentials [10|. 



Since the DKP equation is being increasing 



tivistic spin-0 and spin-1 bosons 
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y used to describe the interactions of rela- 



m 



m 



13, 



la 



19|, it would be interesting 



to probe whether the DKP equation is amenable to exact solutions in the framework of the 
asymptotic iteration method (AIM). This is precisely the aim of this paper. 

In the next section, we explain the AIM briefly and show how to solve a second-order 
homogeneous differential equation. Then, we introduce the DKP oscillator and Coulomb 
problems and obtain their exact eigenvalues and eigenfunctions. In section |3 we present 
the solution of the anharmonic oscillator as a nontrivial example within the perturbation 
theory. Finally, in the last section, we provide our summary and conclusion. 

II. BASIC EQUATIONS OF THE ASYMPTOTIC ITERATION METHOD (AIM) 

We briefly outline the asymptotic iteration method here; the details can be found in 
references 8|, |9|, |lO| . The asymptotic iteration method was proposed to solve second-order 
differential equations of the form 



y" = Xo{x)y' + soix)y 



(1) 



where Ao(x) 7^ and so(x), Ao(x) are in Coo(a,b). The variables, so(x) and Ao(x), are 
sufficiently differentiable. The differential equation ([T]) has a general solution ^ 



y{x) = exp 



adx 



// . , // 



C2 + Ci / exp / Xo{x ) + 2a{x )dx dx 



if n > 0, for sufficiently large n, 



An A„_i 



a 



where 



An — A„ + Sn-l + AoAn-1 



and 



S'ti — Sfi-l + ■^oAn- 



(2) 



(3) 



(4) 



The quantization condition of the method together with equation (^ can also be written 
as follows 



6{x) = \n+l{x)Sn{x) - \nix)Sn+l{x) = 



(5) 



For a given potential, the idea is to convert the relativistic wave equation to the form 
of equation ((T)). Then, sq and Aq are determined and Sn and An parameters are calcu- 
lated. The energy eigenvalues are then obtained by the condition given by equation Q. 
However, the wave functions are determined by using the wave function generator, namely 
exp (— f^ adx ). 

In this study, we seek the exact solution of DKP equation for which the relevant second 
order homogenous linear differential equation takes the following general form. 
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1-6x^+2 X J^ 1-6x^+2- 

If this equation is compared to equation (^, it entails the following expressions 

(m + 1] 



An 



ax 



N+l 



6x^+2 



X 



soix) 



wx 



N 



1 - 6x^+2 



(6) 



(7) 



while the condition Q yields for A^=-l,0,l,2,3,.... 



w'^{-l) = n{2a + 2bm + {n + 1)6) (8) 

w™(0) = 2n (2a + 26m + (2n + 1)6) (9) 

w;^(l) = 3n (2a + 26m + {3n + 1)6) (10) 

w';^{2) = 4n (2a + 26m + (4n + 1)6) (11) 

w™(3) = 5n (2a + 26m + (5n + 1)6) (12) 
... etc 

Hence, these formulae are easily generalized as; 

<W^MiV + 2)-n(n+ P-;;)^^ + ^" ) (13) 

The exact eigenfunctions can be derived from the following generator: 

Using equation (jSj) and equation (fTj), the eigenfunctions are obtained as follows; 

yo{x) = 1 



,.(x) = C.(iV + 2)V(. + l)fl-^^i^±^x--^^ + ^^i^^±^li^±^.^(--^) 



cr a(cr + l) 

a(a + l)(a + 2) 
^^^"^ = -''' (iV + 2)-3 

X ^ _ 3fe(P + 3) ^jv+2 , 3b^(p + 3)(p + 4) ^2(7V+2) _ ft=^(/> + 3)(p + 4)(p + 5) ^3(^^,A 

V ^ c7(a+i) p(p+i)(p+2) ; 

... etc 

Finally, the following general formula for the exact solutions y„(x) is acquired as; 

l/„(x) = (-l)"C2(iV + 2)"(a)„2i^i(-n,p + n;a;6x^+2) (^^^ 



whprp (a) -I>±2) n- - "^^+^+3 j _ (2m+l)b+2a 

wnere [a)n- r(o- ' ^ - Ar+2 ^^'^ P - (Ar+2)fe • 



III. DKP HARMONIC OSCILLATOR 



In this section, the Duffin-Kemmer-Petiau formahsm ll|, |l3| is briefly sketched and the 
DKP oscillator is solved using AIM. Generally, the first order relativistic Duffin-Kemnier- 
Petiau equation for a free spin zero or spin one particle of mass m is 



(c/3.p + mc^)^ = ih0 



,0^ 

dt 



where j3^ (/i= 0, 1, 2, 3) matrices satisfy the commutation relation 



Z?'^/?''/?^ + /3^/?''/?^ = g^'^'p^ + g^'^lS^' 



(16) 



(17) 



which defines the so-called Duffin-Kemmer-Petiau (DKP) algebra. The algebra generated 
by the four (3 matrices has three irreducible representations: a ten dimensional one that 
is related to S=l, a five dimensional one relevant for S=0 (spinless particles) and a one 
dimensional one which is trivial. 

In the spin-0 representation, (3^ are 5x5 matrices defined as (z = 1, 2, 3) 



/3" 



e q\ . / p' 
Ot y V -p^ 

with 0, 0, as 2 X 2, 2 X 3, 3 X 3 zero matrices, respectively, and 



e 




P' 




P' 




P' = 



0-1 







For spin one particles, jS'^ are 10 x 10 matrices given by 



(3' 



^ 0^ 
0^ I I 

0^ I 
0^ 

where s, are the usual 3x3 spin one matrices 



V 



7 



/ 



/3^ = 



Ci 

0^ -iS: 

ej 

0^ -isi 



\ 











ei 



1 



62 



1 



es = 1 



(18) 



(19) 



(20) 



(21) 



I and are the identity and zero matrices, respectively. While the dynamical state ipoKP 
is a five component spinor for spin zero particles, it has ten component spinors for S* = 1 
particles. 

For the external potential introduced with the non-minimal substitution 

p ^ p — imujrfv (22) 



2 



where u is the oscillator frequency and 77'' = 2(3^" - 1, the DKP equation for the system is 

[c/3.(p - imurfr) + mc^] ^ = ihp^^ (23) 

In the spin zero representation, the five component DKP spinor 



^(r) = I with tpupper = I I and tjJlower 






(24) 



so that for stationary states the DKP equation can be written as 

mc^0 = E>^ + ic(p + imujv).K 

mc^(f = E(j) (25) 

mc^A = ■ic(p — imujr)(f) 

where A is the vector {Ai, A2, A3). 
The five-component wavefunction ip is simultaneously an eigenfunction of J^ and J3 

2 I Wupper I _ I -^ Wupper \ ^ T/ T _, i \ | ^"PP^*" 1 /cyf-\ 

Wlower j \ l-^ "I" '-'J Ylower j \ Wlower 

I Y upper \ I J-'^Y upper \ t, r I r upper . /ni-T\ 

Js I I ^ 1 1^1 ' 

Vlower / \ {L3 + S^jyJiower ) \ Wlower 

where the total angular momentum J = L + S which commutes with /9°, is a constant of 
the motion. 



For S = DKP oscillator eigenstates problem, the most general solution for a central 
problem jll| is presented as follows 



where 



i'jMir) 



^ Fnj{r)YjMm ^ 

Gnj{r)YjM{n) 
\^T.LHnjL{r)Yli,{Q) J 



aj 



v/(J+l)/(2J + l), 0=\/^/(2J + l) 



(28) 



(29) 



F„j(r)=F(r), G„, = G(r), H^,j,j±i{r) = H±,{r) (30) 

ipjM of parity (—1)'^ is inserted into equation (|^ and the following equations are found. 



EF = mc^G 



(31) 



^c ( \ -— F = mc^Hi 

dr r hi aj 



(32) 



, , d J muor\ ^ 1 o rr 

dr r h J Qj 



(33) 



d J + 1 mujr \ /■ f ^ "^ mur 

dr r h J \dr r h 



jH_i = ^ {mc^F - EG) (34) 



From the above equations, if equations (jHH) to (jH^ are inserted into equation (0^. the ho- 
mogenous second order differential equation for the DKP harmonic oscillator |11|] is obtained 

as; 



d^ {E^-m^c^) 3muj m^u?r^ J {J + I) 



dr'' 



+ 



(he) 



+ 



h 



ff 



F{r) = 



If we define Eeff 



[E'^-m^c 



m 



2 h ^ and A; = ^, equation (jSSI) becomes 



(35) 



^ + E.,r-er^-'-^^\F(r) 



dr"^ 



(36) 



The asymptotic iteration method requires selecting the wave function as follows 



F{r) = r-^+^e-2>'' f{r 



equating it into equation (J36|l leads to 



+ {E,ff -3k~ 2kJ) f{r) = 



(37) 



(38) 



dr'^ \ r J dr 

where Aq = 2 (fcr — ^^) and Sq = 3k + 2k J — E^jf. By means of equation (@)), we may 
calculate \n{r) and s„(r). This gives: 

2 ( A;r - -+ ^ 



J + 1 



Ac 

Sq = 3k + 2k J — Eeff 



+ 2kJ - E^ff + A\ kr 



J + 1 



Ai = 5k + 2- 



si = 2 {3k + 2kJ - Eeff) I kr - 
Aa = -4 — — + 2{kr- 



4{k + ^^^1 + (3A; - 2k J - Eeff) + 2 ( kr - ^^^ 



S2 = {3k + 2kJ -Eeff) 
etc 



7k + 4^^4^ + 2k J - Eeff + A(kr ^ ^ ^ 



(39) 



Combining these results with the quantization condition given by equation (jSJ yields 

^ (Ee//)o = 3A; + 2kJ (40) 

^ {Eeff)^ = 7k + 2kJ (41) 

■ (Ee//)2 = llA; + 2/cJ (42) 
. . . etc 



Sq _ Si 

Ao Ai 

S^ _ S2 

Ai A2 

S2_ _ S3 
A2 A3 



When the above expressions are generalized, the DKP oscillator eigenvalues turn out as 



{Eeff)n = k{4n + 3 + 2 J) 



(43) 

If one inserts the values of k and Eeff into equation ()43|) . the relativistic energy spectrum 
of DKP oscillator becomes 

1 



2mc'^ 



{E. 



NJ 



m 



,2 4 



c^) = Nhw 



(44) 



where A^ is the principal quantum number defined as A^ = 2n + J. Our result is in agreement 
with the result of reference [11| for the same potential. 

As indicated in Section El we can construct the corresponding eigenfunctions by using 
the wave function generator given by equation (J14p and equation (|39j) where we obtain A 
and s values. Therefore, similar to equation (fT3)). the wave function /n,(r) can be written: 

/„(r) = (-1)" C22^{a)n iFi {-n, a; kr^) (45) 

F{r) ensues right away in the following form: 

F(r) = r'+'e--^'^" [(-1)" C22"(a)„ iF, (-n, a; kr^)] (46) 

where 

^ = ^— and {a)^ = -^^r^ (47) 

Using the wave function F{r), the wave functions G{r), Hi{r) and iJ_i(r) can be easily 
obtained by using equations (jHH) to (jMj) . 

IV. DKP COULOMB POTENTIAL 

We now apply the AIM method to the bound state problem of a spinless charged pion 
(tt^) in the Coulomb field of a nucleus. If we use following ansatz: 

mc'±E V X ^ 2 . , , ^, . 27E 

a± = , 7 = aZ, A^ = — , K = -—\jv[i^& — F^, 4 = — ;^, p = nr (48) 

Aic mc he nnc 

the system of coupled equations for the Coulomb potential becomes 

(dF _J +^-r\ _ ^ rr 

\dp p J k\ 



dp P J \ dp p 



Eliminating Hi and H^i in favor of F, the second-order differential equation for the 
Coulomb potential becomes 

£F{p) fi 1 J(J+l)-72- 



dp^ ' \p A p^ 

Let the radial wave function be factorized as: 



F{p) = 



F{p) = p^+^e-5^/(p) 



where 



Equation (jK^ becomes 



A 



h^^'^l^' 



Y 



(52) 



(53) 



(54) 



d'fip) (p-2A-2)d/(p) (A + 1-0 



dp2 



dp 



P 



f{p) = 



(55) 



which is now amenable to an AIM solution. In order to find the exact energy eigenvalues, 
we define Aq and sq as 



An 



(p - 2A - 2) 



So 



(A + 1-0 



P P 

Using the quantization condition given by equation (j3)), the ^ values take the form 



(56) 



6=A + 1, 6 = A + 2, e3 = A + 3, 



(57) 



which can be generalized as 

^„ = A + n (58) 

Inserting ^ and A in equation ()48|) and defining the principal quantum number as n = n + J, 
we obtain the exact bound state eigen-energies: 



Enj = mc 



1 



{azy 



(n-J-\ + ^(j+iy-{azy 



(59) 



which is in agreement with the results of the references 12, lldl for the same potential. 
The binding energy Bnj can be calculated from Bnj = mc^ — Enj- 
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We can also construct the corresponding eigenfunctions using AIM as 

Up) = {-If C2(a)„ iFi i-n, a- p) (60) 



which gives 



where 



F{p) = p^^'e-'^^ [(-1)" Cicr)^ ,F, {-n, a; p)] (61) 



a = 2A + 2 and (^)„ = ^^^^ (62) 

Other components of the the wave functions (G(p), Hi{p) and H^i{p)) can be obtained 
through equations (PUJ) to (jHUI) using F{p). 

V. ANHARMONIC OSCILLATOR 

In this section, we present the application of the Asymptotic Iteration Method (AIM) to 
non-trivial problems. We have thus chosen a vector potential of type: 

f/v = r2«, ^ = 2,3... (63) 

Taking ^ = 2, the second-order DKP equation becomes as follows 

^.nr)+[ ^, S-^jfW=0 (64) 

In order to solve this equation with AIM, we propose the following wave function to transform 
it to an equation similar to equation (jT}: 

F(r) = e-^'^^''f{r) (65) 

where (3 is an arbitrarily introduced constant to improve the convergence speed of the 
method. We take it /5=5 as in reference J2y] to compare with their non-relativistic re- 
sults for a similar problem. By taking h = c = m = 1 and J = (s-state) for simplicity and 
inserting this wave function into equation ()64|). we obtain 

i^/(r) = {-E^ + 2Er^ + (3 + l- P^r' - r«) /(r) + 2 /? r^/(r) (66) 

11 



which can be now solved by AIM. Here, the so{r) and Ao(r) are as follows 

So(r)= (-^2 + 2 Er^ + /5 + 1-/5 V-r*^), Xo{r) = 2(3r (67) 

In order to obtain the energy eigenvalues from equation (jMj) . using equation Q, we obtain 
the Sfc(r) and Xk{r) in terms of So(r) and Ao(r). Then, using the quantization condition of 
the method given by equation |S1 we obtain the energy eigenvalues. This straightforward 
application of AIM gives us the energy eigenvalues, however, we have observed that the 
energy eigenvalues oscillate and do not converge within a reasonable number of iteration. 
The sequence appears to converge when the number of iterations fc <~ 30, but then it begins 
to oscillate as the iteration number k increases. This result violates the principle behind 
the AIM; as the number of iteration increases, the method should converge and should not 
oscillate. We have noticed that the first reason for the oscillatory behavior is the r* term 
and the second but less serious reason is the E'^ term. 

Therefore, in order to overcome this problem, we have used a perturbation approach 



within the framework of the AIM, similar to reference 



2l|. In order to apply the perturba- 



tion, we introduce a parameter 7 for so(r) in equation (jHH): 

so(r) = {-E^ + 2Er^ + (3 + l + -f{-(3^r^ - r^)) (68) 

7 is an artificially introduced perturbation expansion parameter and at the end of the cal- 
culations, it will be seen that it is equal to 1. After this, equation (0) becomes 

h{x, 7) = Afc+i(a;, 7)5^(0;, 7) - Xkix, 7)sfc+i(x, 7) = (69) 

If we expand 5(x, 7) near 7=0, we obtain the following series: 

6k{x,-f) = 6k{x,0) + - — — +— ^-^ — +— ^-3^ — +... (70) 



1! (97 



=0 3! dj 



7=0 



7=0 2! dy 
According to AIM, the quantization condition 6k{x,'-f) must be equal to zero: 

It is also suitable to expand the energy eigenvalue E, 

En = El + ^El + ^^El + ^^El + ^'Et, + ... (72) 

12 



En expansion terms can be obtained by comparing the terms with the same order of 7 in 
equations (f7T|) and (f?^ . Hence, it is clear that the roots of 61 (x,0)=0 give us the main 
contribution energy terms E^ and the roots of Sj^ (x, 0)=0 give us the first correction E^ 
and so on. 

After we apply this perturbation approach, we have obtained the ground and the first even 
excited state energy eigenvalues. The results are presented in Tables HI and HT] respectively 
for the ground and the first even excited state eigenvalues. In the first column of Table HI 
we present the Eq (unperturbed), the second column Eq which is the first correction and so 
on. We have used the perturbation up to 5*^ term, but one can use higher terms to improve 
the results. However, the effect becomes smaller as it can be seen from tables. In the last 



20| for the same 



column of Table HI we show the non-relativistic results of the Fernandez 
potential to compare with our results. For these calculations, we have observed that the first 
term (E^) in the expansion ()72|1 converges around ^=30 iterations, however, the correction 
terms require higher iterations and start to converge around k=50. 

In Table HTl we show the first even excited state energy eigenvalues. Again, the pertur- 
bation is calculated up to S*'^ and the first term (Eq) converges around k=35 iterations, 
however, the correction terms require higher iterations and we have run them up to k=50 
iteration. 

VI. CONCLUSION 

This paper has presented a different approach, the asymptotic iteration method, to the 
calculation of the non-zero angular momentum solutions of the relativistic DufRn-Kemmer- 
Petiau equation. Exact eigenvalues and eigenfunction for the relativistic Duffin-Kemmer- 
Petiau oscillator and Coulomb problems are derived easily. The advantage of the asymptotic 
iteration method is that it gives the eigenvalues directly by transforming the second-order 
differential equation into a form of y" =Xo{r)y' + SQ{r)y. The exact wave functions are 
easily constructed by iterating the values of sq and Aq. We have also shown how to solve 



13 



the non-trivial problems with the help of the perturbation theory within the framework of 
the asymptotic iteration method. The method presented in this study is general and worth 
extending to the solution of other interaction problems. 
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TABLE I: Ground state energy of the anharmonic oscillator where k is the iteration number (n = 0, 
h = c = m = \ and J = (s-state)) 



k El 


El 


El 


^0^ 


Et 


El Eo Eo [20j 


5 2.478891 


-0.481521 


-0.171317 


-0.087565 


-0.038408 


-0.030214 1.669866 


10 2.477792 


-0.485884 


-0.158642 


-0.080739 


-0.054255 


-0.036055 1.662217 1.325073435 


15 2.477837 


-0.485459 


-0.159187 


-0.082888 


-0.052218 


-0.035973 1.662112 1.147766154 


20 2.477839 


-0.485450 


-0.159249 


-0.083021 


-0.051830 


-0.035991 1.662298 1.072223000 


25 2.477838 


-0.485452 


-0.159247 


-0.082987 


-0.051875 


-0.036052 1.662225 1.062711298 


30 




-0.159246 


-0.082983 


-0.051885 


-0.036069 1.662203 1.060482716 


35 




-0.159246 


-0.082984 


-0.051885 


-0.036062 1.662209 1.060372025 


40 








-0.051884 


-0.036060 1.662212 1.060362059 


45 










-0.036061 1.662211 1.060362077 


50 










1.060362091 


55 










1.060362091 


60 










1.060362090 


65 










» 55 55 


70 


•>1 








)) 5) 55 
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TABLE II: First even excited state energy of the anharmonic oscillator where k is the iteration 
number (n = 2, h = c = m = 1 and J = (s-state)) 

k E2 E2 E2 E2 E2 E2 E2 

5 5.698344 -1.478344 0.027555 -2.574579 -1.656472 3.568314 3.584818 
10 5.370588 -0.911718 -0.040349 -0.433823 -0.040637 -0.068883 3.875178 
15 5.413951 -0.974055 -0.267240 -0.202274 -0.037847 -0.001956 3.930579 
20 5.415995 -0.992837 -0.286535 -0.104029 -0.066451 -0.078990 3.887153 
25 5.415458 -0.991076 -0.277803 -0.121674 -0.072980 -0.054753 3.897172 
30 5.415453 -0.990550 -0.276740 -0.127545 -0.073126 -0.040145 3.907347 
35 5.415460 -0.990604 -0.277153 -0.126693 -0.071566 -0.043842 3.905602 
40 5.415460 -0.990626 -0.277226 -0.126297 -0.071080 -0.046266 3.903965 
45 5.415460 -0.990623 -0.277201 -0.126333 -0.071315 -0.045760 3.904228 
50 5.415460 -0.990622 -0.277194 -0.126358 -0.071409 -0.045369 3.904508 
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